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Toroidal magnetic fields subject to the Tayler instability can transport angular momentum. We show that the Maxwell 
and Reynolds stress of the nonaxisymmetric field pattern depend linearly on the shear in the cylindrical gap geometry. 
Resulting angular momentum transport also scales linear with shear It is directed outwards for astrophysical relevant 
flows and directed inwards for superrotating flows with AQ,/AR > 0. We define an eddy viscosity based on the linear 
relation between shear and angular momentum transport and show that its maximum for given Prandtl and Hartmann 
number depends linear on the magnetic Reynolds number Rm. For Rm ~ 1000 the eddy viscosity is of the size of 30 in 
units of the microscopic value. 
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1 Introduction 

Magnetic field configurations in astrophysical objects occur 
in wide variety. They play an essential role in the forma- 
tion of flow pattern, are source and result of dynamo pro- 
cesses and enhance transport processes, where especially 
transport of angular momentum is of interest here. Cru- 
cial for magnetic field evolution is the interplay between 
toroidal and poloidal fields. In a differentially rotating envi- 
ronment for instance a weak poloidal field is wind-up into a 
strong toroidal one. A strong toroidal field on the other hand 
can become unstable due to the Tayler instability (Tayler 
1957, Tayler 1973, Vandakurov 1972) or the Azimuthal mag- 
netorotational instability (AMRI, see Rudiger et al. 2007). 
As result also new components of a poloidal field are gener- 
ated (Gellert, Rudiger & Fournier 2007). To study in more 
detail the occurrence of the Tayler instability (TI) and es- 
pecially the connected transport of angular momentum, we 
use a model in a cylindrical gap, that starts with an already 
dominating toroidal field and look for the fields instability 
in a differentially rotating environment. 
Tayler instability mainly leads to nonaxisymmetric fields. 
This gives raise to the idea it could be part of a new dynamo 
mechanism caused only by the magnetic field (Spruit 1999). 
Even if there is a ongoing discussion if such a dynamo 
could really work (Braithwaite 2006, Zahn, Brun & Mathis 
2007,Genert, Rudiger & Elstner 2007,Rincon et al. 2008), 
the TI might play an essential role for mixing of chemi- 
cals (Rudiger et al. 2008), transport of angular momentum 
and resulting flattening of the rotation profile of stars. Es- 
timates for angular momentum transport based on the so- 
called 'Tayler-Spruit' dynamo are in use as input for star 



evolution models (Maeder & Meynet 2005, Yoon, Langer & 
Norman 2006). Putting such estimates on more robust ba- 
sis is the aim of this work. It uses nonlinear 3D simulations 
of the Tayler instability in a cylindrical gap to compute the 
actual angular momentum transport 

T^{u'i,u'^-—B'^B'^), (1) 

MOP 

where field fluctuations are defined as nonaxisymmetric 
contributions of the whole velocity and magnetic field in 
cylindrical coordinates (i?, 0, z). Further on it is shown that 
T scales linear with the differential rotation, that is angular 
momentum transport can be expressed in terms of a posi- 
tive turbulent or eddy viscosity i^t- We find that for given 
shear and magnetic field strength always exists a maximum 
value of the eddy viscosity, which grows linearly with the 
magnetic Reynolds number Rm. 

2 Basic equations and numerics 

We use the magnetohydrodynamic Fourier spectral element 
code described by Gellert, Rudiger & Fournier (2007). With 
this approach we solve the 3D MHD equations 

dtu + (it • V)u = -Vp + V^u - 



Ha 

— (rotB)xB, (2) 
Pm 



dtB - ^V'B 
Pm 



divM = 0, 



rot(M X B) 
divB = 



(3) 
(4) 



for an incompressible medium in a Taylor-Couette system 
with inner radius Rin and outer radius i?out- Free parame- 
ters are the Hartmann number 
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(7) 



i) = -Rout/-Rin = 0.5 and fj,n = fiout/i^in. Following the 
Rayleigh stability criterion 9fl(i?^fi)^ > a configuration 
with fiQ > 0.25 is hydrodynamically stable. This is always 
fulfilled in our computations. In addition to the shear flow a 
toroidal external field i?ext = (0,5^,0) is applied. The 
second component has a radial profile similar to the 
flow with a component that is connected to a current only 
within the inner cylinder wall and a second with a current 
also within the fluid domain: 



(8) 



Thus the external field than can be characterized by the 
number 



Ms 



Bout 



asRout + bsRout 



(9) 



-Bill flS^in + bsRin 

The current-free field is given hy fiB — 0.5. Values below 
lead to a field connected with constant negative current jz 
within the flow and larger values to a positive current. The 
Tl occurs in the case that 



(10) 



In the following fiB is fixed to the value //b = 1. This means 
a nearly constant field along radial direction, which is unsta- 
ble after the stability condition (fTOl i. 

The initial magnetic field in the simulations consists of 
small random perturbations (white noise with an amplitude 
of lO^^SP ). 



LU 

m 



1000 



500 



and the magnetic Prandtl number Pm = h'/rj. Here v is 1500 
the viscosity of the fluid and rj its magnetic diffusivity. The 
Reynolds number is defined as Re = r2iiii?iiil?/:/ with D — 
Rout — Rin (unit of length) and the angular velocity of the 
inner cylinder Qi^. Unit of velocity is i^/D and unit of time z 
the viscous time /v. ^ 
The solution is expanded in M Fourier modes in the o 
azimuthal direction. This gives rise to a collection of merid- > 

HI 

ional problems, each of which is solved using a Legendre cc 
spectral element method (see e.g. Deville et al. 2002). Ei- 
ther M = 8 or A/ = 16 Fourier modes are used, three ele- 
ments in radial and eighteen elements in axial direction. The 
polynomial order is varied between iV = 10 and N = 16. 
With a semi-implicit approach consisting of second-order 
backward differentiation formula and third order Adams- 
Bashforth for the nonlinear forcing terms time stepping is 
done with second-order accuracy (Fournier et al. 2004). 

Boundaries at the inner and outer walls are assumed to 
be perfect conducting. In axial direction we use periodic 
boundary conditions to avoid perturbing effects from solid 
end caps. The periodicity in z is set to F = QD. 

The simulations start with an initial angular velocity 
profile 

n{R)=a+^, (6) 
the typical Couette profile, with 



Pm=0.1 
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Fig.l Stability map for Tl with Pm = 0.1, /to = 0.5 
and /iB = 10. There exist two separated regions where the 
instabihty occurs. 

3 Results 

3.1 Tayler instability 

The Tl needs basically two conditions to be fulfilled to 
work. This is on one hand a strong enough magnetic field for 
given values of fi^ and jl^. And, on the other hand, it needs 
a not to high Reynolds number, because otherwise the in- 
stability is suppressed by the shear. In Riidiger et al. (2007) 
stability maps for different Prandtl numbers resulting from a 
linear analysis are presented. The structure of all these maps 
for Pm < 1 is roughly the same. For the case of Pm = 0.1 
such a stability map is shown in Fig. [T|with values result- 
ing from computations with the full 3D code. It differs only 
slightly from the linear results. Remarkable is the separation 
of the unstable region by a stable island. Nevertheless the 
appearance of the Tl is in both regions the same. The most 
unstable mode is the nonaxisymmetric mode m = 1. And in 
both regions exist a upper boundary for the Reynolds num- 
ber where the instability disappears. With smaller Pm the 
critical Reynolds numbers grow. For Pm = 1 it is roughly 
Re = 450 for Hartmann number Ha = 150, for Pm = 0.1 
already Re = 1500. 

3.2 Definition of eddy viscosity 

For the calculation of the Maxwell and Reynolds stress we 
define mean fields as quantities averaged in azimuthal di- 
rection: 

{u) = ^jud4>, {B) = ^j,Bd<l>. (11) 

This regards nonaxisymmetric parts as fluctuations u' = 
u — (it) and B' — B ~ (B) respectively. The normalized 
equation for angular momentum transport (eq[T|| then reads 



r:=(«>;)-^(BiBi). (12) 

The question arising now is whether angular momentum is 
transported outwards, e.g. T > for dil/dR < 0. If this is 
true, an eddy viscosity i/j- can be introduced by defining 

T^-ijtR^ (13) 
dR 
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Fig. 2 Radial profile of angular momentum transport T 
for several differential rotation profiles //ji- Reynolds num- 
ber is Re — 30 and Hartmann number Ha = 190. T is 
mainly positive for subrotation dil/dR < and becomes 
negative for superrotation dQ/dR > 0. 



DC 
DC 



C/5 

o 



o 
4 










Ha^/Pm <B' B' 

H ( 
_ _ _ 


> 

» 



-4 






' ' ' 




















" 1 


> 
I. 


















-8 








T " 



0.25 0.5 0.75 1 1.25 1.5 



ROTATION RATIO 

Fig. 3 Scaling of angular momentum transport. Maxwell 
and Reynolds stress with the shear, expressed as rotation ra- 
tio fin of inner and outer cylinder Positive T means trans- 
port of angular momentum outwards. 



with positive vt- 

Radial profiles of T for Pm = 1, Rc = 30, Ha = 190 
and fiB — 1 found in our nonlinear simulations are shown 
in Fig. 12] We find indeed that the value of T depends on 
/tjj. The profile is everywhere positive for fiQ < 0.65 and 
everywhere negative for fin > 1.1. For values in between 
there are regions where T is positive or negative. Taking 
a radially averaged value of T shows, that this average is 
positive for fin < 0.75 and negative for larger fin- 

It further on reveals a linear dependence between T and 
fin in the whole range of 0.3 < fin < 1-5 (see Fig.|3]l. The 
kinetic part decreases with increasing fin, it is zero for solid 
rotation as one would expect. Contrary the magnetic part in- 
creases with fin and is positive. The linear fit crosses zero 
at the Rayleigh line. There the pure hydrodynamic instabil- 
ity dominates and makes numerical checks impossible. The 
effect, that T is not zero for solid rotation but has a contri- 
bution T* flis called A-effect (Rudiger 1989,Kichatinov 
& Rudiger 1993). Its appearance her is not yet investigated 
in detail. 

Because of the linear scaling of T with fin, the introduction 
of an eddy viscosity as defined by eq. (fTSl l is legitimate. If 
we substitute Rdfl /dR by using fin as 



dn 

dR 



(1 - m) 



Rflu 



(14) 



in eq. (fTSI ) and using the slope of the linear function 
T = T{fin) plotted in Fig. [3] we get an eddy viscosity in 
units of the microscopic viscosity of v^/v = 0.39 for the 
example with Ha = 190 and Re = 30. This is rather low 
and due to the fact that Re is chosen to be very small. In the 
astrophysically more relevant regime, where Re is (much) 
larger than Ha the eddy viscosity increases rapidly too. 
Before describing the scaling of with increasing 
Reynolds number and decreasing Prandtl number a remark 
regarding the size of the Maxwell and Reynolds stress is ap- 
propriate. As shown in Fig. [3] both values are of the same 
order With increasing Re and Ha the magnetic part be- 



comes more and more dominating. So the turbulent viscos- 
ity is mainly due to the magnetic field action. On the other 
hand, the velocity correlation stays on this low level. Argu- 
ing that the mixing of chemicals in stars is connected only 
to the kinetic part of the correlation tensor, mixing due to 
the TI would be rather weak (Rudiger et al. 2008) and in 
agreement with observations of lithium concentration in so- 
lar type stars (Brott et al. 2008). 

3.3 Reynolds and Hartmann number dependence of 
eddy viscosity 

To study the influence of turbulent viscosity for instance in 
star evolution models, one would like to have a relation at 
hand connecting the size of v^/v with given (magnetic) 
Reynolds and Hartmann numbers. This are usually large 
values and are not reachable with our 3D simulations, thus a 
direct computation is impossible at the moment. Neverthe- 
less, in a first step we want to give a scaling relation of 
depending on Re and/or Ha for 0.1 < Pm < 1. 
As already mentioned the TI for given Hartmann number 
occurs only in a certain range for Re. Along a line of con- 
stant Ha starting from the lower critical value of Re the in- 
stability first becomes stronger in a sense, that the angular 
momentum transport increases. It reaches a maximum and 
fades away towards the upper critical Re. This behavior is 
shown in Fig.|4]for Pm — 1 and Ha = 190. The fit with a 
second order polynomial suggests a T ~ Rm^ scaling for 
the increasing part. 

The calculation of i^x from T using eq. ( flj] ) gives the 
values plotted in Fig. |5] The dependence seems to be the 
same as for T: i^t ~ Rm^ if far enough from the up- 
per stability boundary. The maximum value in this example 
(Ha = 190) for the eddy viscosity is max(t^T/'^) = 12.3. It 
grows further when Ha is increased and reaches an absolute 
maximum of the order of 30 for Ha « 500 for Pm = 1. 
For Prandtl numbers below Pm = 1, this tendency to reach 
a saturation value is not yet visible in the with our code 
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Fig. 4 Angular momentum transport for Pm = 1, Ha = 
190 with fiB ~ I and fifi — 0.5. The rotational stabilization 
for large Re leads to the existence of a maximum value. The 
fitted line, a parabola, suggests a relation T ^ Rm^. 
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Fig. 5 Same as in Fig. |4]but for the normalized eddy vis- 
cosity I V. 

achievable range for the Reynolds and Hartmann number 
respectively. The several computed values are listed in tab. 
[T] Because the Reynolds number grid is rather coarse, the 
given maxima might be slight underestimations of the real 
maxima, which reside between the given value and the first 
lower point towards the upper stability boundary. 

The next question we attack is if it is possible to find a 
relation between the eddy viscosity maxima and the param- 
eters they depend on, and how such a relation will look like. 
The first idea is that such a relation should reflect the in- 
fluence from all free parameters, thus a quantity depending 
on magnetic field, shear and Prandtl number. Thereby on 
aspect simpUfies the problem: we deal only with the max- 
ima. For them it is equivalent to use Re or Ha, both are 
not independent. So the relation in question might hide the 
dependence on the third parameter Two well-known param- 
eters, the Lundquist number S and the magnetic Reynolds 
number Rm combine the influence from two physical prop- 
erties of the system. Both are first candidates to check. The 
plot of V'y{S) does not reveal a direct relation between eddy 
viscosity and Lundquist number and is not shown here. 

On the other hand a plot of j^t as function of Rm shows 
indications of the relation in question. Using the values of 
V'yIv listed in tab. [Ugives the data points plotted in Fig.|6l 
Those points reveal a linear dependence 

j/t/j^ = 0.026 Rm + 0.45 (15) 



Table 1 Eddy viscosity and the belonging Reynolds 
numbers for which V'^jv reaches its maximum for given 
Prandtl number and Hartmann number 
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Fig. 6 Linear dependence of eddy viscosity maxima on 
the magnetic Reynolds number for the investigated Prandtl 
number range 0.1 < Pm < 1. 

between eddy viscosity maxima and magnetic Reynolds 
number independent on Prandtl number The one excep- 
tional point much below the line for Rm = 320 is due to 
the fact, that this point, belonging to Ha = 100, is to close 
to the left stability boundary (see fig. [T] for Pm = 1 the 
Hartmann number range is very the same). Lets remind the 
reader that vt here means the maxima of the eddy viscos- 
ity produced by the TI for given Prandtl and Hartmann (or 
equivalent Reynolds) number. As visible in Fig.|5]the linear 
relation i^t (Rm) is not true in general for arbitrary combi- 
nations of Hartmann and Reynolds numbers. 

The derived relation ( fTSl ) agrees also quite well with the 
resulting value for the eddy viscosity described in subsec- 
tion 13.21 using the linear dependence between T and the 
differential rotation parameter /isi- Indeed, going to higher 
Hartmann numbers leads to an increase of the eddy viscos- 
ity to a value of i^t/i^ = 1-1, which is nearly exactly the 
value relation ( fTSl l predicts for Rm = 30. 

The scaling with Rm reflects the raising influence of the 
magnetic field fluctuations for faster rotation. It is derived 
for nearly constant magnetic fields {fiB = 1) and a Prandtl 
number between Pm = 0.1 and Pm = 1. If we assume it 
holds also for smaller Prandtl numbers, one could expect for 
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instance in young, fast rotating stars with Rm « 1000 a sig- 
nificant transport of angular momentum and eddy viscosi- 
ties of the order of 30. For a sun-like star with Rm w 150 it 
would be a weak effect with ut/u = 

4 Conclusions 

Strong toroidal fields under the influence of a differentially 
rotating medium can become unstable to nonaxisymmetric 
perturbations. This instability, the Tayler instability, appears 
only if the shear is not to strong. In cylindrical gap geome- 
try for Prandtl numbers 0.1 < Pm < 1 there is always an 
upper boundary for the Reynolds number above which the 
instability is suppressed. The smoothing action of differen- 
tial rotation dominates. Additionally it needs the crossing 
of a threshold for the strength of the magnetic field, for a 
nearly constant toroidal field with piB = 1 this is of the 
order of Ha « 100. If our astrophysical object in mind pro- 
vides such requirements, the TI leads to nonaxisynmietric 
magnetic field pattern. 

Using the nonaxisynmietric parts as fluctuations, the re- 
sulting angular momentum transport scales linear with the 
differential rotation, it is positive for astrophysical relevant 
configurations with /to < 0.75, e.g. angular momentum is 
transported outwards. This gives raise to the introduction of 
an eddy viscosity which causes the additional turbulence in 
the underlying medium. 

For a fixed magnetic field there is always an upper and 
lower critical Reynolds number in between the instability 
can be observed. And there is also a maximum value of an- 
gular momentum transport and eddy viscosity. Taking the 
maxima of the latter we find its values follow a Unear de- 
pendence on the magnetic Reynolds number. A direct influ- 
ence of the field strength is hidden behind the fact, that with 
the Reynolds number also the Hartmann number is fixed for 
the point where the eddy viscosity has its maxima. 

The linear relation VTiJiva) is derived for a simplified 
model with fixed toroidal field and without any influence 
from density stratifications for rather high Prandtl numbers 
0.1 < Pm < 1. Although, imagine it hold also for Prandtl 
numbers one or two orders of magnitude less, this relation 
could be used to estimate eddy viscosities usable as input for 
star evolution models. For instance for a young, fast rotating 
sun it gives values around vt^/u = 30. 
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